Stability of Resonant Opto-Mechanical Oscillators 
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We theoretically study the frequency stability of an opto-mechanical radio frequency oscillator based on 
resonant interaction of two optical and one mechanical modes of the same optical microcavity. A generalized 
expression for the phase noise of the oscillator is derived using Langevin formalism and compared to the phase 
noise of existing electronic oscillators. 
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I. INTRODUCTION 

The opto-mechanical oscillator (OMO) generates spectrally 
pure radio frequency (RF) signals \Xr^\ due to pondero mo- 
tive interaction between photons and phonons. Phase noise 
and linewidth are the main characteristics detremining the per- 
formance of an oscillator and its practical usefulness. It was 
shown that the OMO linewidth can be small enough to be ul- 
timately described by a Schawlow-Townes-like formula {T\. 
A Leeson model |8] of the phase noise far from the carrier- 
for a radiation pressure driven OMO was presented in [9]. In 
this paper we derive a generalized formula for the OMO phase 
noise that takes into account the noise of the light pumping the 
resonator. Using the formula, we analyze Allan deviation of 
the oscillator frequency as well as its linewidth. The predicted 
performance of the OMO is then compared with the perfor- 
mance of existing electronic oscillators. 

Opto-mechanical oscillation can be described as a strongly 
nondegenerate parametric process in which a pump photon is 
transformed to a photon of lower frequency (Stokes photon) 
and a phonon |10]. This occurs when the pump power ex- 
ceeds a certain threshold determined by the loss of the system, 
and the coupling efficiency between the light and the mechan- 
ical modes. In the case of microcavity-based OMO, the pro- 
cess is allowed if the resonator has at least two optical modes 
with frequency equal to the frequency of a mechanical mode. 
This condition satisfies energy conservation. The momentum 
conservation (phase matching of the parametric interaction) is 
met by the requirement of spatial overlap between the optical 
and the mechanical modes that ensure the convolution integral 
of the electric fields and the mechanical displacement, as well 
as mechanical strain, is nonzero lUlll . 

Generally, the radiation by a laser is assumed to be phase 
insensitive. It means that the phase of the emitted light in 
a laser does not depend on the phase of the pump, and that 
there is no phase matching conditions (preferred directional- 
ity) in the system. An OMO, even one based on Brillouin 
lasing Ill2l[l3ll . has different properties. The mechanical qual- 
ity factor is high enough to endow the generated phonons with 
some well defined phase. The phase of the Stokes light emit- 
ted by the OMO depends on the phonon phase, so the ampli- 
fication of the Stokes wave becomes phase sensitive and the 
phase fluctuations of the optical pump leak to the phase fluctu- 
ations of the Stokes light. In what follows we develop a model 
for an OMO and find the phase noise of the generated signal. 

It was shown that the linewidth of an OMO increases if gen- 
eration of an anti-Stokes optical sideband is allowed by the 
system architecture tTnlOJ. In this paper we analyze an ideal 



case where only two high quality (Q-) factor optical modes 
and a single mechanical mode interact. The optical pump is 
resonant with higher frequency optical mode and the Stokes 
optical sideband is generated in the lower frequency mode. 
The bandwidth of the optical modes is assumed to be much 
smaller than the mechanical frequency. This kind of interac- 
tion has been realized in oscillators based on stimulated Bril- 
louin scattering llalOl and surface acoustic waves 0141 [iSll . 

We write three coupled Langevin equations for the modes 
and solve them analytically to evaluate the phase noise of the 
OMO. To obtain a consistent result, we take into account the 
phase noise of the laser used for pumping the OMO. We show 
that one limitation of an OMO is that the frequency stabil- 
ity of the mechanical mode is determined in the same way as 
the electronic quartz oscillators. The fluctuations of the ther- 
mal bath as well as the thermodynamic noise and drift of the 
mechanical mode limit the frequency stability. The pump- 
ing light is used as the power source in an OMO, similar to 
the electric power source in electronic oscillators and, there- 
fore, it is unlikely that the OMO will outperform its electronic 
analogs, unless the frequency bandwidth of the optical modes 
is significantly smaller than the frequency bandwidth of the 
mechanical mode. An advantage of an OMO, on the other 
hand, is in its potential long term stabilization via stability 
transfer from the optical frequency domain. 

The paper is organized as follows. The equations are pre- 
sented in Section II, and their solutions are described in Sec- 
tion III. The OMO phase noise and linewidth are analyzed and 
compared to similar parameters of existing quartz oscillators 
in Section IV. 



II. BASIC EQUATIONS 

The triply-resonant opto-mechanical interaction is de- 
scribed by equations 



A = -TaA - igCB + Fa, 
B = -TBB-igC^A + FB, 
C = -rcC-igB^A + Fc. 



(1) 
(2) 
(3) 



where A, B, and C are the slowly-varying amplitudes the 
pump (optical), the Stokes (optical mode red shifted with re- 
spect to the pump), and the signal (mechanical) fields; Ta, 
Tb, and Tc are the linear resonant terms of optical and me- 



chanical modes respectively 

r^ = i{uJa - Wo) + 7 + 7cQ, 

Tb = i(wh - w_) + 7 + 7c6, 

Tc = i{uJc-UJM) +1M, 

7 and 7m are the intrinsic decay rates of the optical and me- 
chanical modes, ^ca and ^cb are optical loading (coupling) 
rates (the loading of the optical modes can be different be- 
cause modes belong to different families); g is the opto- 
mechanical coupling constant. 



and 



5 = Wo 



K,h 



2m*L'^ojc 



(4) 



Kf^ is the correction coefficient showing that radiation pres- 
sure results not only in a change in the size of the resonator, 
but also in its index of refraction through strain, m* is an ef- 
fective mass of the mechanical mode, L is an effective spatial 
parameter of the mode. 

The terms Fa and Fb represent Langevin forces with two 
uncorrelated parts arising from the internal and coupling loss 
of the modes 



FA^F,A + FrA, {FA)^e''^^A 



2Pjc. 



(5) 



{F,A{t)FUt')) 

{FrAit)FUt')) 



hojo 

27ca<5(t-t'), (6) 

2j5(t-t'), (7) 



FB^F,B+FrB, {Fb)=0, (8) 

{FcB{t)F^^s{t')) - 2j,bS{t-t'), (9) 

{FrBit)F^^sit')) = 2j6it~t'), (10) 

where (...) stands for ensemble averaging, and P is the 
power of the external optical pump of the mode A. 

The Langevin force describing the thermal fluctuations of 
the mechanical system is defined in the similar way 



{Fc) - 0, 
{FcWFlit')) 



2jM{nth + l)5{t-t') 



(11) 
(12) 



where fith = [exp{hujc/kBT) — 1]^^ is the averaged number 
of thermal phonons leaking from the thermal bath of tempera- 
ture T to the mechanical mode, fc^ is the Boltzmann constant. 
We present the slowly-varying complex amplitudes as 



B= \B\e"^^, 
C= \C\e">'^, 



(13) 
(14) 
(15) 



and derive two sets of equations for the amplitude and phase 
parts of the complex amplitudes from Eqs. (I2]|3} 



l^l = -(7 + 7ca)|A|-.g|C||B|sin0 4 

\{FA)\cOii{(j)FA - </'a) + Fat, 

\B\ = -(7 + 7c6)|B| +g|C||A|sin0 4 
\C\^-lM\C\+g\B\\A\sm^ + Fcr; 



F, 



Br; 



(16) 

(17) 
(18) 



\C\\R\ 

(j)A^ ~{uJa-UJo) - g COS(f)+ (19) 

\{Fa)\ ■ ,, I \ . ^Ai 

——— sm{(j>F^ ~^'^'~^Ta\' 

(/)j3 = -(u;b-cj-) -g COS0+-— |-, (20) 



1^1 , , 

; , ^ \B\\A\ Fc^ ,-., 

(pc = -{^c-i^Mj- g ,„, cos<l)+j^; (21) 



\B\ 
Fa 

\C\ 



where 



T>A- (PB - (PC, 
1 

2* 



Fat ^ ^{FcAe~"^^ + Flj^e'-t 



FBr ^\{FcBe-'^^ + Fl^e'*^) + 
^{F,Be-'^-+Fl^e'^-), 
Fcr = \{Fce-'^^+Fle'^^), 
Fm = ^{F,Ae-'^^ - Flj^e'^^) + 



FB^ = ^iF.Be-'*^ - Fl^e'^^) 



1 



-{FrBe-'^" ~ Flge'^-), 



1 
2i 

Fc^ = ^iFce-'^^-Fle'*^). 



(22) 
(23) 

(24) 

(25) 
(26) 

(27) 

(28) 



Equations ( I16H21| | completely describe the triply-resonant 
opto-mechanical interaction. 



III. SOLUTION 



Using the set of equations for the phase ( fT9ll2Tl i we get the 
equation for phase difference (p 



'{oJa - UJb -Wc) +.gcos( 
- sin{(pFA - (f>y 



B\\A\ ^ \C\\A\ \C\\B\ 



\C\ ' \B 

\{Fa)\ ... : ._^FA^ Fb^ Fa 

^^sm(0^,-0^) + ^-^-^, 



1^1 



ip = ~{uja ~LOb- LOc) -.9 sin?/; 



\C\ 



\B\ 



\A\ 



, \{Fa)\ . ,. . .^FA^ Fb, Fc, 

+ ^^sm(0^,^0^) + — - — -—. 

Next, we introduce an expectation time-independent value for 
the phase difference, (tp), and fluctuational time-dependent 



(29) 



We introduce %p ^ (p — Tr/2 and rewrite Eq. (l29t as 

\B\\A\ , \C\\A\ \C\\B\ 



(30) 



part, S4>, so that -0 — {"4^) + 5(i>, and get 



5sin(V') 
I(^a)| 



|B||A| , |q|A| |C||i?| 



sm 



1^ 

5sin(V') 
{Fa)\ 



\C\ \B\ 

\B\\A\ , \C\\A\ 



\A\ 



\C\\B\ 



(31) 



\C\ 



\B\\A\ 
\C\ 



\B\ \A\ 

, \C\\A\ \C\\B\ 



\A\ 



^Fa 



\B\ 
Fm 

\A\ 



1^1 
Fbi Fa 



(32) 



where 5[. . .] means a deviation from the expectation value. 

Using Eqs. ( I16H21I ) we find several relationships for the ex- 
pectation values of the oscillator parameters in the steady state 



\Bl 
\A?-^ 

UJb - 



IM 



1 + lcb 


' 


IM 


r^P 


1 + lcb 


g' 


^- 7 + Icb 



- ujm 



IM 



a^W 



(33) 

(34) 

(35) 
(36) 



The expectation value of the amplitude of the field in 
the pumped mode increases below the oscillation thresh- 
old with increase of the pump power, in accordance with 
\A\ = \Fa\/\Ta\, and then stays constant, in accordance with 
Eq. dll. 

For the sake of simplicity, we assume that the Stokes 
sideband has much lower power compared with the pump, 
A| ^ I i?|, and that the system is triply -resonant (the expecta- 
tion values of the frequencies are resonant with corresponding 
modes). We note that 



\B\\A\ , Idl^l \C\\B\ 



\C\ 



\B\ 



IM 



1 + lcb 



\A\ 

1 + lcb 
IM 



IFbI 



\B\ 



(37) 



IM 



ir, 



1 + lcb 9 



so that for the case of relatively weak Stokes sideband and 
all-resonant tuning 



9 



\B\\A\ \C\\A\ \C\\B\ 



\C\ \B\ 

Using the assumptions and 

l^ll^l 



1^1 



IM +1 + Icb- (38) 



5 

\{Fa 



-gsin(V') 



1^1 



•^Fa 



\C\ 



we find 

\C\\A\ 
\B\ 



\C\\B\ 



\A\ 



(39) 







\{Fa)\ 
\A\ ^ 



where we took into account that, in accordance with Eq. (|3Tt 
and Eq. ( |36] |, {4>Fa ~ 'Pa) =0 for the resonant tuning. 



We obtain a set of linear equations consisting of Eqs. (fT9] l 
and(|32]i 



H + iiM + 1 + icb)54' = 

(50^^_J^^)K£p+^A. Fs. Fa 



\A\ \A\ \B\ Id' 

bA = [dcpFA - ocpA) 1^1 + -r^r- 



(40) 



(41) 



The condition of small Stokes sideband can be written in the 
form 



\{Fa) \ 
\A\ 



= l + lc 



Substituting Eq. ( |42] | into (|40l i and ( 1411 ) we derive 

FA^ 



^A + (7 + 7. 



ca)"VA 



\AV 



(42) 



(43) 



bB + S(t>c + ilM + 1 + lcb){HB + 50c) = (44) 
Fsi , Fa 



ilAI +1 + lcb)5<j)A + 



\B\ \C\ 



The third equation for phase deviations can be derived from 
Eqs. (|2i) and (EB: 



..W-^-50c/^ + ^" 



7 + 7c6 



IM 



1 + lcb \B\ V IM \C\ 



IM Fbi jl + lcbFc- 



(45) 



We solve the set of linear differential equations (1431145b us- 
ing Fourier transform, e.g. 



FB^ 

Fa 



fBii^)e~ 



/cj(a;)e" 



fduj 
"2^' 
fduj 
'2^' 



(46) 
(47) 



where fBi{^) and fci{oj) are the Fourier components of the 
noise. 



{FB^{t)FB^{t')) = 1 (7 + 7,,)5(t - t') , 



(48) 



{Fa{t)Fa{t')) = 1m (nth + ^j S{t - t'), (49) 

{fB^ii^)fB^{L0')) - ^i^ + lcb)S{uJ + Uj') , (50) 

{fa{^)faW)) = 27r7M Uth + ^ S{io + uj'). (51) 

Readout of the OMO signal is accomplished by tracking 
thradio frequency beat note produced by the optical pump and 
the generated sideband on a fast photodiode, so the phase 
and frequency of the measured oscillation signal are deter- 
mined by the argument of the product of optical amplitudes 
(AB*). The phase noise of the signal is given by difference 
6(j)B " S(j)A- On the other hand, the mechanical frequency can 



be read using electronics means, e.g. a capacitive displace- 
ment sensor Neglecting by the electronics back action, we 
can estimate the phase noise of the signal evaluating d(j)c'- We 
find expressions for Fourier amplitudes of these parameters 



S4>Ai^) 



fA 



1 






(52) 
(53) 



(5(/)a(w) 



-iuj + 7Af + 7 + 7cfc 

-iuj + 7m /b£ 

-iuj{-iuj + 7M + 7 + 7cb) \B\ 

7 + 7cb fch 

'iuj{-iuj + 7a/ + 7 + Icb) \C\ ' 

IM 



5(t)c{uj) 



-ioj + 7Af + 7 + 7c6 
7a/ fsi 



6a (^) - (54) 



'iuj{-iuj + 7a/ + 7 + Icb) \B\ 

-ioj + 7 + 7cb fct 

'iuj{-iuj + 7A/ + 7 + 7cb) \C\ 



Equations ( |53] | and ( l54l l can be used to find the phase noise 
of the OMO. For example, single-sideband phase noise Cd^) 
of the mechanical oscillation is defined as, 

oo 

(50c(i)50c(i-T))- /■ /:^(,,)e-*^. (55) 



Using definition 



Ki) = 



/OO 7 

<50c(^)e-^"*-^, (56) 



we find 



(57) 



1m 



7^ 
lea 



w^ + (7m + 7 + 7cb)^ ^2 + (7 + 7ca)2 

7 + 7cb , 



*-'i' 



ill 



w2[a;2 + (7A/+7 + 7cb)2] 2|i?|2 

^^ + (7 + 7cb)^ 7a/ /- 1' 

oj^[uj^ + (7A/ + 7 + 7cb)2] |C|2 1^"*" + 2; ' 

where £i„ stands for Fourier frequency dependent input phase 
noise of the pump laser. This value is usually much larger than 
the contribution from the quantum white noise of the laser 
Using similar reasoning we obtain 



^n-h — 



(58) 



w^ + ill 



7^ 
lea 



OJ^ + (7m + 7 + 7cb)^ ^2 + (7 + 7co)2 

7 + 7cb , 



r 



^' + ih 



a;2[a;2 + (7A/ + 7 + 7cb)2] 2\B\^ 

il + lcb)'^ IM f. 1 

L^2[^2 + (^^,^ + ^ + ^^^)2] |C|2 1^"*'^ + 2 



It is also useful to write an expression for the phase noise of 
the optical Stokes mode, determining the stability and spectral 
purity of the Brillouin laser 1I12I1 



Cb 



(59) 



(7 + 7cb)' 



7^ 

lea 



w^ + (7m + 7 + 7cb)^ ^2 + (7 + 7cq)' 



.r 



^"+lM 



7 + 7cb 



^^[i^^ + ilM+l + lcbY] 2|i?|2 

(7 + 7cb)^ 7m_ f- 

c^2[a;2 + (7A/ + 7 + 7cb)2]|C|2 P'' 



Formulas ( |57] |. (ISST i and ( |59] l can be simplified further using 
the ratio between the photon number in the Stokes mode and 
the phonon number in the mechanical mode ( [33] l, and the ex- 
pression connecting the number of pump and Stokes photons 
and output power of the pump and Stokes light 



\A\^ = 



27c 



P 



(7+7ca)^ ?iwo' 



\B\ 



Pbo 



27cb?iajo 



(60) 



where we assumed that the carrier frequency of the Stokes 
light is approximately equal to the frequency of the pump 
light. Finally we get 



i^a-b 



(61) 



^' + ill lla 


r 


^2 + (7A/ + 7 + 7cb)^ a;2 + (7 + 7^^ 


i) 


a;2 + 272^(^4,, + 1) 7cb(7 + 7cb) 


huJo 


[w2 + (7A/ + 7 + Icbf] ^2 


Psout 


Cc = 




ill lla 


-rC,,, 



(62) 



^^ + {iM + 7 + 7cb)^ a;2 + (7 + 7^^)2 

ill {^th + 1 + wV[2(7 + 7cb)^]} 27cb(7 + 7cb) fu^o 



Cy,= 



[^2 + (7A/ + 7 + 7ch)^ 

(7 + 7cb)' 



w^ 



Psout 

(63) 



^ + (7m + 7 + 7cb)^ a;2 + (7 + 7^0)2 
uj"^ + 2jljinth + l) Icbil + lcb) ^iwo 



+ {iM + 7 + 7cb)^] tu^ 



Pbo 



(64) 

Equation (I6TI 1 gives a complete description of the phase 
noise characteristic of the radio frequency photonic oscillator 
based on demodulation of the light output of the OMO on a 
fast photodiode. Equation (|62| | shows the limiting phase noise 
of the oscillating mechanical mode that could observed, if the 
mechanical oscillation is measured using an external devise 
that does not disturb the system. Equation (|63] | describes the 
phase noise of the SBS laser These expressions can be used 
to find the linewidth and Allan deviation of the correspond- 
ing OMO signals. By definition, the Allan variance of the 
frequency of the oscillator is given by 



a\T) = 



4w2£sin4(tjT-/2)dw 
~l4 (wr/2)2 2^' 



(65) 



Let us consider the case where the OMO signal is retrieved 
by demodulation of the light leaving the resonator on a pho- 
todiode. The single sided power density of the phase noise is 
related to the linewidth of the oscillator Aiy as C{uj — s> 0) = 
l-K^vjuP'. We find for the cases of relatively low and rela- 
tively high-Q of the mechanical mode 

i ca 



\1 + 7caj 



lch{l + lch)(nth + 1) 






Bout 



Al^a-b\-tM<.l 



ill 



Ic 



-Aiy„ 



7 + 7cb 



7mK'i + 1) 



(7 + 7c.)M7 + 7c.)^ ^^™^ 

hujQ 



t^Pbou 



(66) 



(67) 



Therefore, in the case of low-Q mechanical mode, the 
linewidth of the radio frequency beat note generated by the 
OMO is determined by the linewidth of the pumping light; 
while in the case of high-Q mechanical mode, the linewidth 
is described by Schawlow-Townes-like formula [TJ and the 
phase noise of the pumping light is suppressed. 



we also show the phase noise of a commercially available 100 
MHz oven controlled quartz oscillator (hne (5) in Fig. ([TJ). 
Apparently, the electronic oscillator has much lower phase 
noise than the OMO, even though thermal drifts resulting in 
the flicker noise are not taken into account in our analysis. 




10° 10' 10' 10' 10" lO'^ 10° 10' 10" 
Frequency f, Hz 



IV. DISCUSSION 

Let us estimate the phase noise of the oscillator and com- 
pare it to the phase noise of an electronic oscillator We as- 
sume that Qm = 10^, lom = 27r X 10^ rad/s, 7^1/ = 27r x 
i4 ^^A /c ., . _ ^^^ = 27r X IQS rad/s. 



500 rad/s, 7 = 27r x 10^ rad/s, 7^;, 

^''punip 



Wo = Svr x 10^'' rad/s, Avpump = 1 kHz, Pboui = 100 /iW, 



and F = 1 mW. 

The light escaping the cavity is demodulated on a photo- 
diode to produce the radio frequency signal. The photodiode 
introduces thermal noise and white shot noise 



CpD = 



2qRpPD + kBT 

Prf 



(68) 



in addition to the phase noise (|6T1 ) coming from the opto- 
mechanical oscillator Here p is the resistance of the photo- 
diode, R — rjq/TiLUo is the responsivity of the photodiode, 
rj is quantum efficiency of the photodiode, q is the electron 
charge, and Pd is the total optical power reaching the photo- 
diode. We assume that Pd = P. The thermal noise depends 
on the ambient temperature T and the power of the radio fre- 
quency signal leaving the photodiode Pjip — 2pR^PPBout- 
To find its value we assume that the resistance at the output of 
the photodiode is p = 50 Ohm, the responsivity of the pho- 
todiode is 0.8 A/W, and the temperature is T = 300 K. For 
these parameters, the expected power of the radio frequency 
signal escaping the photodiode is Ppp — 6 /iW. 

Adding Cpo and Ca-b we plot the spectrum of the single 
sideband phase noise (line (1) in Fig. ([U). It is dominated by 
the phase noise of the pump laser. If a very narrow linewidth 
laser is used instead, the phase noise of the OMO would be 
given by the convolution of curves (4) and (3) (phase diffu- 
sion of the oscillator as well as shot noise). In the same figure 



FIG. 1: Phase noise of the opto-mechanical oscillator (1) character- 
ized with parameters defined in the text. Curve (2) defines the con- 
tribution from the phase noise of the pumping laser. Curve (3) stands 
for the thermal and white shot noise defined in Eq. ( 1681 . Curve (4) 
shows phase diffusion of the opto-mechanical oscillator. Curve (5) 
stands for a commercially available 100 MHz ovenized quartz oscil- 
lator. 

Let us evaluate the sources of frequency drift of the opto- 
mechanical oscillator We note that, according to (l35l l and 
energy conservation law, loq = a;_ + ujm. 



1 + lcb 



ujQ-uj- = COM = (69) 

IM 



7 + 7c& + 7a/ 



7 + 7c6 + 7m 



(cjo - Wh). 



The OMO frequency drifts if the frequencies of the optical, 
ujb, and mechanical, tUc, modes drift. 

It is possible to lock the frequency of the pump light to the 
frequency of the corresponding optical mode, ojq ~ oJa- Then 
the OMO frequency becomes dependent on the parameters of 
the microcavity only: uja — (^b as well as w^- This is similar 
to the case of a quartz oscillator, except the OMO frequency 
depends on both the eigenfrequency of the mechanical and 
optical modes. The dependence of the oscillation on the drift 
of the mechanical mode can be suppressed if the bandwidth 
of the optical mode is much smaller than the bandwidth of the 
mechanical mode. 

A clear advantage of the OMO over an electronic oscillator 
is that it can be readily stabilized by locking the pump fre- 
quency to the frequency of the coiTesponding cavity mode, 
ujQ — oja, and then locking the temperature of the cavity to 
a thermally insensitive optical reference line, for example, an 
atomic transition. Such a stabilization of the temperature will 
stabilize the long term drift of the oscillation [|l6ll . 



V. CONCLUSION 

We have theoretically studied the phase noise of a triply- 
resonant opto-mechanical oscillator based on a nonlinear op- 
tical microcavity. The oscillator generates Stokes optical pho- 
tons and mechanical phonons out of photons of a coherent 
pumping light. We have shown that the spectral purity of the 
opto-mechanical signal is primarily limited by the phase noise 
of the pump laser. The overall short term performance of the 



oscillator is expected to be worse than the performance of a 
conventional electronic oscillator. An important advantage of 
the opto-mechanical oscillator is the possibility to optically 
stabilize it. With a proper electronic locking scheme it is pos- 
sible to transfer frequency stability from the optical domain to 
radio frequency domain. If stabilized to an atomic transition, 
an opto-mechanical oscillator can be made to have better long 
term stability compared to an oven controlled quartz oscilla- 
tor 
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